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Abstract 



One loop QCD corrections to hadronic W decay are calculated for arbitrary W po- 
larizations . The results are applied to W pair production and decay at LEP200. 
We focus on the corrections to angular distributions with particular emphasis on 
azimuthal distributions and correlations. The relevance of our results to the exper- 
imental determination of possible nonstandard triple gauge bosons interactions is 
discussed. 



1. Introduction 



The LEP200 experiment at CERN will start data taking above the WW threshold 
in 1997. While much of the emphasis of the LEP200 program has been on the 
possible discovery of the Higgs and the precise determination of m\y, there remains 
the important aspect of probing the triple gauge boson couplings which are an 
essential feature of the non-abelian gauge symmetry of the standard model. Of 
course, there is already indirect evidence for the existence of triple gauge boson 
couplings from the high precision of the data at LEP1 []IJ|, but extracting bounds on 
these couplings requires unavoidable model dependent assumptions. Even though 
the statistical accuracy of the LEP200 data will not be as impressive as at LEP1, it 
will be sufficient to remove any doubts concerning the structure of the pure gauge 
sector in the standard model. 

It is expected that at threshold, radiative corrections may play a significant role in 
W pair production. Electroweak corrections have been known for some time and 
there has been some progress towards including finite W width effects, which arise, 
for example, from the interference between the processes e + e" — > /1/2/3/4 with and 
without intermediate Ws ||. In this paper we present the one loop perturbative 
QCD corrections in the case where at least one of the Ws decays hadronically, 
W — > qq'(g), which to the best of our knowledge have never been calculated along 
the lines we follow ||J]. 

In principle it is possible that higher order corrections modify the standard model 
predictions in such a way that they mimic the existence of nonstandard triple gauge 
boson couplings, if not taken into account. Therefore, our results should be com- 
pared to possible effects from nonstandard couplings as parameterised in j5| |J. It 
is conceivable that nonstandard effects enter at the few percent level and are thus 
of the same order of magnitude as the QCD effects discussed in this paper. 

It would be possible to completely neglect higher order QCD corrections by studying 
only final states in which both the Ws decay leptonically, however much information 
concerning the triple gauge boson vertices will be inevitably washed out due to 
incomplete kinematical reconstruction, quite apart from large statistical errors due 
to low event rates in this channel. Therefore hadronic decays are important. For 
definiteness, we shall mostly consider in this article the process e + e~ — ► W + W~ — > 
l + vqq'{g) where only the W~ decays hadronically, although our result can easily be 
extended to e + e~ — ► W + W~ — > q^q-^q^- Our treatment is along the lines of |3| 



modifying the discussion of differential distributions sensitive to non standard triple 
gauge boson vertices to include the effects of higher order QCD corrections. Let us 
stress once again that the W's are assumed on-shell throughout the paper (as in 
H) so that our result does not take care of finite width effects whatsoever. 

In section 2 we shall discuss "inclusive jet angular distributions" (i.e. the matrix 
element for W — > qq'g is integrated in such a way that it adds to the lowest order 
cross section e + e~ — > W + W" — > l + uqq') , but for completeness, we include the 
fully differential results in the appendix. In eq. [25| we present a nice compact 
formula, which summarises the most important piece of our result. In section 3 we 
shall discuss in detail some phenomenological consequences of our result on various 
differential distributions. 



2. A Detailed Description of the Calculation 



Non expert readers might be tempted to think that QCD corrections to W — > qq' 
and consequently also for e + e~ — > W + W~ — > l + uqq', are given by a correction 
factor (1 + cHg/it). However, this is true only as far as the total rate is concerned, 
but is not the case if one considers some of the angular disributions relevant for 
probing triple gauge boson vertices. The corrections to the angular distributions are 
related to QCD corrections to hadronic Z decay M, but are more complicated due 
to the larger number of particles in the final state. Correspondingly, more angles 
are necessary to describe the cross section. 

Following the treatment of || the cross-section can be written as 
da(e + e- -> W + W~ -> l + vm') 



d cos dd cos 6id<ftid cos 



E Faba>b>(s, cos ??)£%(#, <P)D A , B ,(7r - B h & 

A,B,A',B',\ 



where the sum runs over the polarizations A,B=L,+,- of the W~ and A',B'=L,+,- 
of the W + and the electron helicity A. D° AB are the "decay functions" of the decay 
W~ — > qq', normalized such that 

+ 1 r 2n AA 4 

dcosd %(^) = -^ (2) 
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They depend on two angles (polar angle 9 and azimuthal angle <fi) which specify 
the direction of the outgoing quark q in the W~ rest frame with respect to the 



direction defined by the W~ motion in the lab frame. The definition of 9[ and (pi 
is similar. The functions F ABA , B , can be constructed from the helicity amplitudes 
for e + e _ — ► W + W~ . They depend on the total e + e~ energy y/s and the production 
angle d of the W's. Since we are interested in QCD corrections to W-decay we shall 
not discuss the functions F\ BA!B , any further. Explicit expressions can be found in 
ref. §. 

For the sake of convenience we reproduce from || the lowest order decay functions, 

D° LL (e,<j>) = sin 2 9 



0L(M = ^(i±«**) 2 

D + _{e,d>) = lsm 2 6e 2i ' l> ^ 



2 

D° ±L (9, 0) = (± cos 9 sin 9 - sin i 
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These results can be derived by contracting the "hadron tensor" H fiy (W~ —>■ qq') 
with the corresponding W polarization vectors e A , Dab = H^e^e*^ . All the rest of 
the decay functions can be obtained by using the relation 

D AB (9,<P)=D BA (9^) (4) 

which is true beyond the leading order. By summing up all diagonal decay functions 
one gets 

Diui = 2 (5) 

Alternatively, this can be obtained from D tota i = J2a b H uu {—g uu + — i V J ^ ) , where 
W denotes the 4-momentum of the W~ . It should be stressed that throughout this 
work, even when discussing the phenomenology in section 3, 9 and are defined in 
the rest-frame of the decaying W and not in the lab-frame. 

The D functions can be decomposed into the sum of a symmetric and an antisym- 
metric part under the simultaneous exchange 9 <-> 7r — 9 and <fi «-> + ir . When 
we discuss the phenomenology of hadronic W decays in section 3 we will consider 
only the symmetric pieces of the D functions, because in hadronic W-decays, quark, 
antiquark and gluon jets cannot be distinguished and therefore the antisymmetric 
parts of the D functions drop out. In this section, however, we shall present com- 
plete results both for the symmetric and the antisymmetric terms of the W decay 
functions. 



Our aim is to calculate the one loop QCD corrections to the decay functions Dab- If 
gluon emission is taken into account, these functions depend on three more variables 
X, Xi and x 2 in addition to 9 and <ft, all defined in the W~ rest frame, as follows: 
x\ and x 2 are the rescaled energies of the quark q (with 4-momentum pi) and the 
antiquark q' (with 4-momentum p 2 ). They can be given in terms of invariant dot 
products as 

2p 2 p 3 = m 2 w {l - x x ) 2p x p z = m 2 w {l - x 2 ) 2p x p 2 = m 2 w {\ - x 3 ) (6) 

X3 is the rescaled energy of the gluon, with momentum p% and with X3 = 2 — x 2 — x\ 
from energy conservation. 9, <fi and x are the angles needed to fix the spatial 
orientation of the triangle which is formed by the vectors p\, p 2 and p% in the W~ 
rest frame (in lowest order there are only two angles 9 and necessary to fix the 
quark(=-antiquark) direction in the W~ rest frame). After integration over the 
additional variables Xi x i and x 2 , the form of the cross section eq. [I] will be left 
unchanged, with 

D AB (9, 0) = D%(9, 0) + C F ^-D\ B (9, 0) (7) 

replacing D° AB in eq. D l AB is the sum of all one loop QCD corrections (virtual, 
soft and hard gluons). We shall actually integrate in such a way that the triangle 
is 'reduced' to the thrust direction, where the thrust is defined in the ordinary way 
except that it is determined in the W rest frame. For example, for three partonic jets, 
q ,q' and g, the thrust direction is just the direction of the most energetic parton. 
This definition is useful, because experimentally quark, antiquark and gluon jets 
cannot be distinguished. The transformation to the W rest frame should be no 
problem, too, because the hadronic part of the event in general can be completely 
reconstructed. 

The calculation of the D\ B (9,<ft) will be described in the following. More precisely, 
we shall calculate the correction to the ratio ® AB , because in this ratio all ultravi- 

Dtotal 

olet, infrared and collinear singularities present in the higher order matrix element 
drop out. This is a consequence of the universality of those singularities, which has 
been known to be true in QCD for many years. Otherwise such singularities would 
appear in intermediate steps of the calculation. For example, one has 

Dtotai(qq'g) = C F —- X l±^l D o (8) 

totally) * 2ir (1 - xx)(l - x 2 ) total K ' 

with singularities for x\ )2 — > 1. Those ratios have the further advantage that the 
virtual gluon exchange corrections drop out (for the case of massless quarks which 



we assume throughout). Writing 



Dtotai = D° total + C F ^-Dl otal = 2(1 + ^) (9) 
Zir it 

(defined including virtual gluon exchange) we have 

^ = ^(l + **") (10) 

Utotal 2 lix ID\ B 

with 

X AB : = 2D\ B - D\ B D] otal 

= g /' dx x f 1 dx.H^qq'g)^*" - D%{-g, v + ^)) 

JO Z7T Jo A-xi 77% (11) 

What is the form of the polarization vectors e^? In the actual calculation we 
worked in the W~ rest system and we choose the quantization axis to be given by 
the direction of the quark momentum, q 1 = 22 tFxi(1, 0, 0, 1) (or, more generally, 
the thrust direction). Furthermore, we defined the coordinate system such that 
q2 = ^2^Xi(l, 0, sin0 12 , cos$i 2 ) is the antiquark momentum, where cos0 12 is given 

by 1 H — — simply from energy-momentum conservation. Therefore we 

had to produce the polarization vectors from the ordinary vectors (0, 0, 0, —1) and 
^(0, ±1, i, 0) by an arbitrary Euler rotation given in terms of 0, 9 and x, 

e L — (0, — sin 9 cos x, — sin 9 sin x, — cos 9) (12) 

e ±i<t> 

e± = — /=(0, — icosx T cos sin x, — isinx ± cos 9 cos x, ±sin0) (13) 
V2 

Note that the angles x an d Q enter the angular correlation structure of the cross 
section for Z — > qqg as well. The generic form of the lowest order LEP1 cross section 
is an expansion in powers of cos 9 

Me+e ] - f ~* m) = ^ (1 + cos 2 9) + -al sin 2 9 + -a% cos 9 (14) 
dcos0 8 4^ 

where the last term is the parity violating vector-axialvector interference term, 

which cannot be measured, because quarks and antiquarks cannot be distinguished. 

For a three jet process one has an expansion 

2nd(T(e + e- -> Z -> ggg) 3 dau(x u x 2 ) 2 n\ , 3 dcr L(xi, x 2 ) . 2n 

; tt; — ; — ; = ; — ; (1 + cos ") + t — ; — ; sm ^ 

dcos9dxdxidx 2 8 dx\dx 2 4 dx\dx 2 

3da P (x 1 ,x 2 ) 3 do-AT^i,^) . 0/1 

+ I j ,7 cos0- ^ sin 20 cos x 

4 dx\dx 2 2y2 dx\dx 2 

3da T (x 1 ,x 2 ) . 2 3 rfo- M (xi,x 2 ) , . 

+ j , , cos 2x sin - sin0 cos x (15) 

4 dx\dx 2 2v2 dx\dx 2 



with ^d^dxf 1 = ^ ^Ixidx^ ( as l° n & as one neglects higher orders, i.e. two gluon 
emission). If one integrates over x, one recovers the structure of the lowest order 



eq. |14|. In that case, for normalized distributions, there are only two independent 



quantities which fix the cross section completely, namely aL = — and 



dtotal dtotal & total 



These quantities, after integration, correspond to L and P eqs. 18 and 21, and have 



been introduced in a different context in the last two references of |IIJ. In hadronic 
Z decays the parity violating piece cannot be determined, because quark and 

G total 

antiquark jets cannot be distinguished. 

In our case we shall find matrix elements which depend on all three angles 0, \ 
and 9 and have to integrate over \ to get the D\ B , eq. |7]. The appearance of an 
additional angle is due to the presence of an additional plane, the plane spanned by 
the momenta of e and W ± . 

In the difference Xab all singularities and contributions from virtual gluon exchange 
drop out. For example, for A=B=L one has 

D LL sin 2 9 a s r 2w dx f , f , ( n Xi + x 2 — 1 



D 



total 
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+ cos 2 x sin 2 ff) — (2— — -| -) sin 9\ 2 sin y sin 9 cos 9 \ (16) 

Xi (l-xi)(l-x 2 ) J 

in case that the thrust (=quantization) axis is given by p\, i.e. X\ > 0:2,3. Note 
that t-. — sm w 1 12 — 7 is integrable. Analogous results are obtained in the other two cases, 

(l-a;i)(l-X2) to to ' 

x 2 > #1,3 and x 3 > Xi£- Integrating these results over the appropriate phase space 
regions M one obtains 

Dll Sin2d <l-3LC F ^) + LC F ^ (17) 



Dtotal 2 2n 2n 

where 

L = 0.4875 (18) 

is the numerical result of the integration over x% and x 2 - Similarly, for the other 
decay functions: 

D+ + + D - = l±^(l - 3LC F £) + 2LC F ^ (19) 
Dtotal 2 2-n 2-n 

D ++ - D__ _ a 



- cos 9(1 + PC F -^) (20) 



Dtotal 27T 



where 



P = -1.340 (21) 

is the characteristic correction for a parity violating contribution as can be derived 
from eq. 36 of the appendix. This and various details concerning the matrix ele- 
ments will be discussed in the appendix. For the phenomenological applications we 
have in mind here, it is not relevant, because quark and antiquark jets cannot be 
distinguished and therefore the parity violating contribution cannot be determined 
experimentally. Note that, in addition, there are other contributions from the ma- 
trix elements, which disappear when the integration over x is performed (see the 
appendix). 



Next, we find 



!±W^(l-3Z^) (22) 

Dtotal 4 Z7T 



± sin # cos #(1-3LC F ^)- sin 0(1 + PCV^) (23) 
2v2 I 2vt 2vt J 



Dtotal 2\[2 I 

Note that we know D to tai from eq. |9|. In any case, the ratios are most inter- 

esting, because they enter in normalized distributions, which are more sensitive to 
anomalous triple gauge boson couplings than the total cross section. 



Forgetting about the parity violating contributions, our results eqs. |I7] - 123] can be 
brought to a very compact form by defining 

bAB = l + ^l-3LcVf (24) 



Using this, eqs. 17-23, can be rewritten as 



A ry 

Dab = D° AB + 25 AB LC F -^ + 0(a 2 s ) (25) 

This is the central result of this article. It means that the QCD corrections can be 
organized such that apart from an overall normalization factor only the diagonal 

decay functions -Dll, D ++ and D are modified. They are modified by a constant 

term 2LCp^. The absolute numerical value of 2LCp a ^ = 0.024 is relatively 
small, but may be relevant in regions, in which the lowest order terms vanish, like 
D\ L at 9 = 0. To see the significance of the term 2LCp^ one may also consider 
D ++ . The constant correction may be written as ~ cos 2 9 + sin 2 9. Comparing with 



the lowest order decay function D° ++ , eq. ||, we see that a purely transverse state 
appears to accquire a longitudinal component, thus changing the relative propor- 
tion of longitudinal and transverse polarisations in the final state. Such effects are 
among the hallmarks of non-standard triple gauge boson couplings. Although the 
absolute numerical value of 2LCp a ^ w is relatively small, it is certainly compara- 
ble in size with possible non-standard effects, which may be even smaller. Concrete 
phenomenological applications will be presented in the next section. 



3. Phenomenological Applications 

Although the QCD corrections are remarkably simple, they modify the shapes of 
all the distributions relevant for the analysis of the triple gauge boson vertex, like 
azimuthal or polar angle distributions of one or both W's. We will present both 
analytical and numerical results for these distributions. Our basic input parameters 
will be \/s = 190 GeV for the total beam energy and a s = 0.11 for the magnitude 
of the strong coupling constant. The latter value is reasonable in view of the LEP1 
result for a s and the uncertainty in the renormalization scale. (The argument of 
a s should be somewhere between mw and y/s.) We shall present all results in the 
form of a ratio of higher order to lowest order prediction, because this emphasizes 
the QCD effects and reduces side effects like the uncertainty in m w etc. 

Let us first discuss the case, in which only one of the two W decays is considered, i.e. 
e + e~ — > W + W~ — > W + qq'. The azimuthal differential distribution may be written 
as § 

da(e + e~ -> W+j-X) _ 3 da(e + e~ -> W+W~) y D 
d cos fid cos 6L d</>_ 87r dcos'd 

where pa,b is the spin density matrix of the W~ as given, for example, in ||. 
j- denotes the highest energetic jet from the W- decay used to define the thrust 
direction. The angles 9 and <fi defining the thrust direction have been renamed to 6L 
and 0_ to stress that they refer to the W~ decay. Note that for hadronic W-decays 
it is usually possible to transform to the W rest system, because the W-momentum 
can be fully reconstructed from the momenta of the decay products. Using our result 



eq. ^5] one can rewrite eq. |26] as 

da( e +e- -> W + gg) = 3 ^(e+e" -> W+W~) a, _ ^ 
d cos $<i cos 6_d(j)- 8tt dcosfi n F 2n 

x 



J2 PAB (D° AB + 2LC F ^6 AB )} (27) 



Carrying out the polar angle integrations one obtains 

da(e + e- ^ W+j^X) 4 A T ^ a s 4,„ n , _ . n , , 

— i - h 4LC F — + -(Re r + _ cos20_ - Im r + _ sm20_) 

ct<7>_ 3 2i 3 

(28) 

for the shape of the azimuthal distribution. Here r + _ is the weighted average of 
P+-, 



J dcosfi 



da 
d cos i? 



(29) 



da | ^ 

In fig. 1 the ratio d ^~ , is shown using the standard model values for r + _ (Im 
r + _ = and Re r + _ = —0.53). Nonstandard triple gauge boson vertex interactions 
modify the 0_ (as well as other) distributions, by modifying r and p, respectively. 
For example, if there are CP violating gauge boson couplings, r and p may acquire 
non- vanishing imaginary pieces, apart from modifications of the standard model pre- 
dictions for the real pieces induced by all nonstandard effects whether CP violating 
or not. 

In fig. 1, $ and 6L have been averaged out completely. Interesting information 
may also be gained by retaining some of the dependence in the polar angles. For 
example, one may consider the <fi dependence of the forward-backward asymmetry 
in in order to keep track of the parity violating terms (~ cosd) in the triple gauge 
boson interaction. The ratio of this distribution for ho and lo is also included in fig. 
1. We conclude that the QCD effects must be taken into account, if one wants to 
check the standard model triple gauge boson vertex to an accuracy of 1%. 

Furthermore, one may examine distributions, in which 0_ instead of 9- is integrated 
out, 



da(e + e- -> W + j-X) 3 da(e + e~ -> W+W 



d cos fid cos 9- I d cos $ 

r 1 

x 



(p++ + P-) 



(1(1 + cos 2 d _) + 2LC F ^\ + PLL {sin 2 + 2LC F ^\ 1(1 + -)(1 - 3LC F '- , 

v I Z7T ' v 111 ' ) 71 ll\ 

(30) 




da | ^ 

Figure 1: The ratio , as a function of 0_. The dependence on 6_ and $ have 
been averaged out (full curve). To obtain the dashed curved we have integrated 
over cosi? in an antisymmetric way, i.e. the integrand is counted negative, if cost? 
is negative. 



In this equation the parity violating piece ~ cos6L has been left out, because it 
cannot be seen in the hadronic W~ decay. Since the relative size of QCD corrections 
is not the same for the transverse and longitudinal part of eq. 30 the shape of the 
6- distribution is modified by the QCD term. In fig. 2 the ratio of the ho to lo 
distribution is shown as a function of both for symmetric and antisymmetric 
integration over Again, the quantitative results contained in this figure will allow 
to discriminate QCD from nonstandard effects. Clearly, the variation of the curves 
in fig. 2 must be taken into account, if one wants to check the standard model triple 
gauge boson vertex to an accuracy of 1%. Note that in all the figs. 1, 2 and 3 below, 
the average of the curves correspond to the overall QCD correction factor 1 + 

If the decay of the W + is taken into account and both W's decay hadronically, 
the situation becomes more complicated. For example, the shape of the double 
differential azimuthal distribution is of the generic form 
da{e+e~ -> W+W~ -> j+j-X) v 

d(f> + d(f>- ABA'B' 

OL OL 

+ 2LC F -^ J2 r A >B>D A , B , + 2LC F -^ ]T r AB D° AB (31) 

A'B' Z7T AB 

where taba'b 1 -, tab and ta<b' are defined in analogy with r + _, eq. |29|. The angles 
0_ and + refer to the decay of the W~ and W + , respectively. If only one of the 




Figure 2: The ratio of the ho to lo as a function of both for symmetric and 
antisymmetric integration over 



W's decays hadronically one obtains a similar result for d ^ ° d( ^ , with only one term 
~ 2LCp... instead of two such terms. Note that in eq. |l]0_ and <p + were denoted by 
4> and (pi, respectively. We present the numerical results on the effect of the higher 
order corrections on the double differential azimuthal distribution , , da , , in table 2 
(for the case that one W decays leptonically) and table 3 (for the case that both 
W's decay hadronically). The entries in the table are the ratio of higher order to 
leading order cross-sections in the respective bins. All angles are in units of 

7r. We see that the variation of the numbers in table 2 is en gross about the same as 
the variations of the curves in figs. 1 and 2, and conclude that in all the cases it is 
necessary to take these QCD corrections into account, if one wants to pin down the 
standard model triple gauge boson vertex to an accuracy of 1-2%. Of course, the 
overall QCD K factor 1 + ^ was known before, but in figs. 1, 2 and 3 and tables 2 
and 3 we see that it is also important to know the deviations from this average. 



Conversely, if one leaves the polar angle dependence and integrates over the az- 



imuthal angles, one obtains 

da(e + e- -> j+j-X) ,3 . 2 da(e + e- -> W+W" 



d cos $d cos 9_d cos 9 , 4 d cos $ 



x + cos 2 + 2LC F ^}{i(l + cos 2 6 + ) + 2LC F ^} 
+ p LLLL {sin 2 9^ + 2LC F ^}{sin 2 £ + + 2LC F ^} 
+ (p LL++ + p LL __){sin 2 + 2LCVg}{^(l + cos 2 9+) + 2LC F ^} 

+ (P++LL + P-LL){sin 2 £ + + 2LCV^}{^(1 + cos 2 0_) + 2LCfc|£}} 

x(l + ^) 2 (l-3LCV^) 2 (32) 

where Paba'b'( s 5 cos ^) denotes the two-particle joint density matrix for pro- 
duction as defined in ref. 0. The complicated looking formulae eqs. 31 and 32 can 
be verified easily by combining the lowest order formulae of reference || with the 
compact form of our higher order result eq. E5L Terms ~ cos 9 + and ~ cos Q_ have 



been left out, because it has been assumed that both W's decay hadronically and 
those terms cannot be detected in that case. 

This case of both W's decaying hadronically poses additional QCD problems, which 
we will not discuss, but which have recently become the focus of some interest || 
||. These have to do with finite width effects, which are partially nonperturbative 
and thus cannot be calculated from first principles, with important consequences for 
the precision to which the Standard Model can be tested at LEP200. 



The distributions eq. 32 may also be considered in case that only one W (e.g. the 
W~) decays hadronically. In that case in eq. 32 one has to add the parity violating 
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Figure 3: The ratio of the ho to the lo of the antisymmetric piece in cos 9 + (coeffi- 
cient of cos# + in eq. 33) as a function of both for symmetric (solid curve) and 
antisymmetric (dashed curve) integration in 



pieces ~ cos# + , because they can be determined from the direction of the / + . 

da(e+e- -> l + vj-X) .3 , 2 da(e + e~ -> W + W~) f, 

~(t)— 1 a J -Up ++++ +P++—+P—+++P- 



d cos fid cos d-d cos 6 '+ 4 d cos $ 

x + cos 2 Q_) + 2LC F ^) hi + cos 2 0+) 
L 2 2tt j 2 

+ Pllll{sui 2 9_ + 2LC F ^} sin 2 # + 
+ (pll++ + Pll — ) (sin 2 9- + 2LCV^}i(l + cos 2 9 + ) 
+ (P++LL + P-ll) sin 2 + cos 2 0_) + ^LC F ^} 

- (P++++ ~ P++— + P—++ ~ P ){^(1 + cos 2 6-) + 2LC F ^}cos6 A 

- {Pll++ ~ Pll — )(sin 2 9_ + 2LC F ^-} costf+j 



x (i + ^)(i_3LC F ^) (33) 

Note that now the factors depending on 9 + do not get the QCD term. In fig. 3 
the ratio of the ho to the lo of the antisymmetric piece in cos#+ (coefficient of 
cos9 + in the last equation) is shown as a function of 9-, both for symmetric and 
antisymmetric integration in The variation of the curves is somewhat larger than 
in figs. 1 and 2, so that essentially the same conclusions can be drawn. 



Conclusions 



In this article we have presented a complete calculation of one loop QCD corrections 
to W-pair production and decay at LEP200, including W-polarization effects. We 
have put special emphasis on the effects on those angular distributions sensitive to 
the structure of triple gauge boson vertices. Our results are of the order of a few 
percent and thus are certainly as important as finite W-width effects and must be 
taken into account for a precision study of the pure gauge sector of the standard 
model. 

This work is part of a larger project, which one of the authors has pursued over the 
last years, namely to calculate QCD corrections to differential distributions using 
the known QCD effects on total rates plus the real gluon matrix elements . This 
method has been applied successfully even to two loop problems |L0| at LEP1. It 
could certainly be used to extend eq. |25| to two loops. We have not attempted this, 
because it is unlikey to be relevant given the magnitude of statistical and systematic 
errors expected at LEP200. 
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Appendix: The complete W decay functions to 0(a s ) 

Before integration over x, x\ and 22 one has a cross section 

da(e + e~ -> W+W~ -> l + vm'g) 
d cos fid cos 9id<pid cos Qd$dx\dx-2,dx 



containing decay functions Dab(9, (f>, X-> %i, x 2 ) which depend on x, xi and X2 in 
addition to and <fi. These will be given in the following. First, for A=B=L, one 
has 

Dll{9-, (j), x, X\, X2) sin 2 9 



-5(1 - x x )5{l - x 2 )5( X ) + C F ^\l(x 1 ,x 2 )(l - 2 sin 2 9 



Dtotal 2 27T L 

+ cos 2 xsin 2 9) + m(xi, x 2 ) sin x sin 9 cos 9 \ (34) 



where D to tai = 2(1 + — ). The functions l(xi,x 2 ) and m(x 1 ,x 2 ) are shown in table 
1. The integral of I over dx\dx 2 is the number L=0.4875 as given in the main text 



eq. |18[ The integral M of m over dx\dx 2 is given in table 1. m does not contribute 



to -5 5^ — i^rr , — ojt because the coefficient of m vanishes when the integral over v 

a cos ttd cos &i dq>i a cos adcp ° /v 

is performed. It contributes only to the azimuthal (=x) dependence in 3-jet decays. 
Secondly, one finds 

(D ++ + D„)(9,(f>,x,x 1 ,x 2 ) l + cos 2 

= = o(l - xi)o(l - x 2 )6(x) 



D 



total 



— Cp^- (m(xi, x 2 ) sin x sin 9 cos 9 + l(xi, x 2 )(l — 2 sin 2 + cos 2 x sin 2 0) j (35) 



= — cos0O(l — £1)0(1 — X2)o(x) 



total 



— Cp^-^p(xi, x 2 ) cos 9 — n(xi,x 2 ) sin x sin j (36) 

The function p(xi,x 2 ) together with its integral P over x\ and x 2 is given in table 
1. It gives a parity violating contribution which is not measurable in hadronic W 
decays. The function n(xi, x 2 ) and its integral N are given in table 1, too. They do 
not contribute to the inclusive distributions discussed in the main text because their 
coefficient vanishes when the integral over x is performed. Furthermore they arise 
from the parity violating part of the W decay and are not measurable in hadronic 
W decays. Finally, one finds 

D + ^(9,(f),x,x 1 ,x 2 ) e 2 ^sin 2 



Dtotal 

2i<f> 



-6(l-x 1 )6(l-x 2 )S(x) 



e 9 a s ( 

+ -ttCft^-I 1{xi,x 2 )(1 - 2cos 2 x + sin 2 cos 2 x 

2 Z7T I 

2 sin 2 + 2i sin x cos x cos0) — m(x\, x 2 ) sin0(z cos x + sin x cos0) \ (37) 



D± L {9,(ft,x,x 1 ,x 2 ) e ±i<j> . 

= — /=(±cos0sm0 — sm0)d(l — xi)d(l — x 2 )d(x) 

D total 2y2 

£ ± i( P Q, s f 

+ —j=CF-^-\^{x\,x-i) sin #(2 cos x sin X ± cos#(cos % — 2)) 

1 1 i 

— x 2 ) sin# + m(xi, x 2 )(± sin x( c °s 2 9 — -) — - cos % cos 9) 



% 1 1 

+ n(xi, X2)(±- cosx — - sin x cos#) > (38) 



Note that D_ + {9, (ft, x, £i, ^2) and Dl^±{9, (ft, x, %i, x 2 ) can be obtained via the rela- 
tion D AB {9, (ft, x, x 1 ,x 2 ) = D* AB (9, (ft, , x, x-y, x 2 ). 
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l(x t , x 2 ) 
n(x 1 ,x 2 ) 
p(xi,x 2 ) 
m(x 1 ,x 2 ) 



q T 

r xl 

Li ) (l-a:i)(l-a: 2 ) l-zi 



_2£2. 

1 1+3:2 l^XQ. 

(l—xi)(l—x 2 ) 1—xi Xl 



-2 

1 



xi 



1+21 2 x\ 

{\—Xi){\—X 2 ) 1-X 2 X 2 



sin^sln & t — 2+212. _ _i 1_ -jl 

lA y (1-X 1 )(1-X 2 ) i-xi 1-X 2 J 



2 £2+1-2/3:3 

\—x\ 



2l 2:3(1—2:1) 



+ 1 



(1 — 21)(1— 3:2) 1— £1 1—2 



integral 
L=0.4875 
N=3.236 
P=-1.340 
M=2.650 



Table 1: Coefficients of the QCD helicity cross sections. The results depend on 
whether quark, antiquark or gluon are used to define the quantization (=thmst) 
axis. 



l<f> + \\<f>-^ 


-0.4 


0.4 -0.8 


0.8 -1.2 


1.2 -1.6 


1.6 -2.0 


-0.4 


1.039 


1.029 


1.038 


1.033 


1.031 


0.4-0.8 


1.037 


1.031 


1.034 


1.035 


1.030 


0.8 -1.2 


1.034 


1.033 


1.034 


1.033 


1.034 


1.2 -1.6 


1.030 


1.035 


1.034 


1.031 


1.037 


1.6 -2.0 


1.031 


1.032 


1.038 


1.029 


1.039 



Table 2: Ratio of ho to lo for the double differential azimuthal distribution , , , 
for various <p + and 0_ bins. <p + and 0_ are measured in units of ir. It is assumed 
here that the W~ decay is hadronic and the W + decay is leptonic. <p + = <pi refers 
to the direction of the charged lepton in the W + decay. 



10+110-^ 


-0.2 


0.2 -0.4 


0.4-0.6 


0.6 -0.8 


0.8 -1.0 


-0.2 


1.076 


1.069 


1.063 


1.064 


1.072 


0.2 -0.4 


1.069 


1.068 


1.063 


1.061 


1.064 


0.4-0.6 


1.063 


1.063 


1.063 


1.063 


1.063 


0.6 -0.8 


1.064 


1.061 


1.063 


1.068 


1.069 


0.8 -1.0 


1.072 


1.064 


1.063 


1.069 


1.076 



Table 3: Ratio of ho to lo for the double differential azimuthal distribution, eq. 31, 
for various <p + and 0_ bins. <p± are the azimuthal angles defined in the decay of 
the W ± and are measured in units of n. It is assumed here that both W's decay 
hadronically. Therefore the average here is (1 + ^) 2 1 + 2^ in contrast to table 
2, where the average is 1 + 24 . 



